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Within the self-consistent constraint Skyrme-Hartree-Fock+BCS model 
(SHF+BCS), we found equilibrium toroidal nuclear density distributions 
in the region of super- heavy elements. For nuclei with a sufficient oblate 
deformation {Q20 < -200 b), it becomes energetically favourable to change 
the genus of nuclear surface from to 1, i.e., to switch the shape from a 
biconcave disc to a torus. The energy of the toroidal (genus=l) SHF-hBCS 
solution relative to the compact (genus=0) ground state energy is strongly 
dependent both on the atomic number Z and the mass number A. We 
discuss the region of Z and A where the toroidal SHF+BCS total energy 
begins to be a global minimum. 

PACS numbers: 2L60.Jz, 25.70.Pq, 25.70.Jj, 27.90.+b 



1. Introduction 

The term doughnut nuclei was coined by J. A. Wheeler [H [2] many 
years ago. In 1970's and 80's the idea of a toroidal geometry in nuclear 
physics was investigated in a framework of the liquid drop model (LDM) 
and shell corrections O H El El E]. In 1990's [El El [TOl [IH [El El and 
more recently [Ml [T5| the toroidal breakup configurations were found in the 
simulations of heavy-ion collisions using the Boltzmann-Nordheim-Vlasov 
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Fig. 1. a) Potential energy curves for toroidal (genus=l) nuclei of different atomic 
numbers Z and mass numbers A as a function of the aspect ratio R/d^ where R is 
the major radius and d the minor radius of the toroid. b) The same as a function 
of the quadrupole moment Q20 defined in Eq. ([T]). 



(BNV) or the Boltzmann-Uehling-Uhlenbeck (BUU) kinetic transport mod- 
els. The stability of toroidal nuclei against a change of the quadrupole mo- 
ment was studied recently in the Hartree-Fock-Bogolibov (HFB) theory [16j 
with the Gogny DIS force p/Tj and in the semiclassical extended Thomas- 
Fermi (ETF) method [18] with the Skyrme SkM* force [19j. 

The purpose of this work is to study the stability of toroidal super-heavy 
nuclei against collective deformations within the self-consistent constraint 
Skyrme-Hartree-Fock+BCS framework (SHF+BCS) with the Skyrme SkM* 
force [19j. Since the classical LDM and semiclassical ETF approaches do 
not include shell effects, we would like to examine quantitatively the role of 
the shell effects in the description of toroidal nuclei. 

The HF+BCS equations were solved using the code HFODD [2(7, that 
uses the basis expansion method in a three-dimensional Cartesian deformed 
HO basis. The details of our model are the same as in Ref. [22j. 



2. Results 

In topology, the genus of a closed surface is equal to the number of 
handles on the surface. For instance, a torus is a surface with genus=l 
while a sphere or disc has genus=0. The shape of a toroidal nucleus can be 
characterized by an aspect ratio R/d, where the major radius R is measured 
from the center of the toroid to the center of the circular meridian and the 
minor radius d is the radius of the meridian. Assuming axial symmetry and 
a uniform mass distribution we can calculate a quadrupole moment Q20 for 
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a given aspect ratio R/d as 



Q20 = ^ I d'rg{fyY2o{cos^) (1) 
A 02/ 2 



4 V Stt cosh 770 



sinh^ ^0 — 5 cosh^ 770 



where in the toroidal coordinates [4J g{r) = yQ{v ~ Vo)i cosh 770 = R/d and 



o3 



Rq is a spherical radius of the same volume V = (47r/3) Rq = 27v'^Rd'^. 

Figure [1] shows the potential energy curves for toroidal nuclei with 
114<Z<138 calculated within the LDM [HIS]. The potential energy curves 
are plotted as a function of R/d and Q20 in the panel a) and b), respec- 
tively. As one can see all toroidal super-heavy nuclei (genus=l) show an 
oblate deformation with (^20 ^ -200 b. In addition, as the atomic number 
Z exceeds 138, the genus=l equilibrium lies at an energy even lower than 
that for the spherical shape. 

In Fig. [2] we display the total binding energy (E^^^) calculated within 
SHF+BSC approach as a function of the quadrupole moment Q20 for ^^^122, 
340^30^ 352^34 364^33 gupcr-hcavy nuclei. The two branches of E^""^ 
corresponding to a compact (genus=0) and a toroidal (genus=l) SHF+BCS 
solutions have been found. Both solutions coexist in the vicinity of Q20 — 
-200 b. But with the further increase of the oblate deformation it becomes 
energetically favourable to change the genus of nuclear surface from to 1, 
i.e., to switch the shape from a biconcave disc to a torus. For Q20 < -250 b 
only genus=l self-consistent solutions exist. 

We also show the self-consistent density distributions obtained for both 
branches with the compact and toroidal topology. Since the SHF+BCS 
model is not restricted to the axial symmetry, the toroidal density dis- 
tributions show the characteristic sausage deformations [4J. The sausage 
instability is responsible for a multifragmentation of toroidal nuclei and in 
the case of 31^122 and ^^^130 the toroidal solutions cease to exist at Q20 — 
-450 b. 

Similarly to the LDM results showed in Fig.dl the genus=l and genus=0 
total energy minima become closer in energy with increasing atomic number 
Z and mass number A. For Z=138 and A=364 the toroidal equilibrium again 
begins to be the global minimum. 

The neutron Fermi energies (A^) oscillate around -6 MeV with a change 
of quadrupole deformation for all nuclei from Fig. [2l The proton Fermi ener- 
gies (A^) stronger depend on Q2o- For the prolate deformations A^ oscillate 
between -2 and -1 MeV but in the region of oblate deformation A^ decreases 
to the value of A"^ at Q20 ^ -400 b. Thus, the toroidal SHF+BCS solutions 
are more stable against /3-decay than the compact (genus=0) solutions. 
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Fig. 2. The total binding energy (E^^^) of the toroidal (genus=l) and compact 
(genus=0) SHF+BCS solutions as a function of the quadrupole moment Q20 for 
340^30^ 352;l34 364;L3g rj.^^ self-consistent density distributions with the 
compact and toroidal topologies are also shown, where a z-axis lies on the surface 
of the page for a prolate deformation and is perpendicular to the page in the case 
of an oblate deformation. 



In conclusion, it appears that our self-consistent SHF+BCS calculations 
for the toroidal nuclei are consistent with the classical LDM and semi- 
classical ETF [18J models as well as with the HFB theory 

The toroidal nuclear density distribution is not a surprising phenomenon, 
but a regular characteristic of the strongly oblate deformed heavy and super- 
heavy nuclei. The doughnut shaped nuclei are not the exception but rather 
the norm in this region of deformations! 
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